Downloaded from rspa.royalsocietypublishing.org 


z 

cc 


2<zui 

W>- LU — 

^ u 

S Q. o3 LO 



Ss 

w (J 

Xo 
H ^ 


PROCEEDINGS THE ROYAL A 
-of- SOCIETY rkI 


MATHEMATICAL, 
PHYSICAL 
& ENGINEERING 
SCIENCES 


On a Method of Finding the Conductivity for Heat 

C. Niven 

Proc. R. Soc. Lond. A 1905 76, 34-48 
doi: 10.1098/rspa.l 905.0003 


Email alerting service 


Receive free email alerts when new articles cite this article - sign up in the box 
at the top right-hand corner of the article or click here 


C/5 

o 

z 

5 

uu 

LU 

U 

O 

O' 

Cl 


u. 

0 


_r o 

< Z 
— CC 

b u* 

5>uj — 
v u 

«Eo.o3i/> 


H 

uj 


-j 

£ 

2 
w u 
Xo 

H t/5 

C/5 

o 


u. 

0 


u 

O 

c* 

Ql. 


To subscribe to Proc. P. Soc. Lond. A go to: http://rspa.royalsocietypublishing.org/subscriptions 


This journal is © 1905 The Royal Society 














Downloaded from rspa.royalsocietypublishing.org 


o 

z 

cc 


2<Zw 


ZQ. 08 (O 


34 Prof. C. Niven. On a Method of [Dec. 5, 





w u 

EO 

H 


C/5 

O 



u 

o 

0 * 

CL 


_T C 3 
< z 

y 5 

I<Zm 

r* z qj 
5 >: uj — 
2 i v u 
zo.o3 io 



w U 


lo 


H ^ 



0^ a Method of Finding the Conductivity for Heat . 

By Professor C. Niven, F.B.S. 

(Beceived December 5, 1904,—Bead January 26, 1905.) 

Introduction. 

This paper contains an account of a method for finding the conductivity for heat, 
especially in bad conductors. The substance is supposed given in the form of a 
cylinder, or rather of two half-cylinders pressed together and heated by a current 
passing through a wire along the axis. When the steady state has been reached, 
the conductivity is given in terms of the difference of temperatures at known 
distances from the axis, and the heat supplied by the wire. The latter can be 
found when we know the current passing through it and the difference of potential 
at its ends. 

The first part of the paper contains a detailed account of the methods employed 
for finding the difference of temperature, and a description of the apparatus used. 
The results of some experiments made with it are also given, and compared 
with those found by other observers. 

The second part of the paper contains a solution of the mathematical problem 
of the diffusion of heat in an infinite solid from a line at which it is supplied at a 
constant rate, and the solution of some other allied questions. 

One result of the investigation suggests a method of finding the diffusivity 
directly, when the substance is of sufficiently great extent. 

Part I. 

1. If heat be supplied uniformly and continuously to a solid within a 
given space at the rate of H units per second, and if a surface S be drawn 
round this space, then when the steady state has been attained and the 
temperature at each point is constant in time, 

H = 

J an 

the integral being taken over the closed surface, h being the conductivity, 
and clO/dn the normal variation of 0 outwards. In particular, if the solid 
he in the form of a circular cylinder, and the heat he supplied by an electric 
current passing through a uniform straight wire coinciding with its axis, 
the equation takes the form 

the enclosing surface being a coaxial cylinder, and Hj the heat supplied per 
centimetre per second. 
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1904.] Finding the Conductivity for Heat . 35 

If we integrate this equation, we find for the difference of temperatures 
#i—at distances a , 6, from the axis, the expression 

Hi log* bja = 2tt — 0 2 ) A. 

Upon this equation we may base a method of finding h, since the heat 
evolved in the wire per unit length is readily found when we know the 
current i (in amperes) passing through it and its resistance; or what comes 
to the same thing, the current and the difference of potentials "V^—Vg 
(in volts) between two points at a known distance apart. These data can 
be got by using, say, a Weston millivoltmeter provided with suitable shunts 
and resistances, and may be obtained with considerable accuracy if the 
instrument be properly calibrated. 

Reducing electrical to heat units the value of the conductivity is given by 

in _ V 3 )log e 5/a 

b 4.2 x2t r(0 x —0g)Z’ 

l being the length of the heating wire, at the ends of which the potentials 
are taken. 

The heating wires used were of platinoid. At first wires of No. 40 
S.W.G-. were used carrying currents of about one ampere, but these were 
afterwards replaced by thicker wires of No. 34, as they were found to give 
too great a heating effect, and caused the wood in which they were placed to 
shrink at the centre, and besides, in addition, were liable to snap when drawn 
tight. The thicker wires gave more satisfactory results. 

2. The temperatures 0 V 0 3 , or rather 9^6^ were measured either by the 
change of resistance of platinum wires embedded in the body parallel to 
the axis, or by thermo-electric junctions formed by soldering together thin 
iron and German silver wires placed in the same positions as the platinum 
wires. These fine wires were stretched round two brass pins so as to be 
at distances about 1 cm. and 3 cm. from the heating wire. 

When solids were tested they were made in the form of two half 
cylinders, one of which was placed in a semi-cylindrical wooden shell, which 
was moved up till the two wires rested on its flat face. If platinum wires 
were used their resistances were now compared. The current is then passed 
through the heating wire, which is drawn tight round a brass pin and fixed 
to lie parallel to the other two wires. 

The second half cylinder and the other semi-cylindrical wooden shell 
were then placed upon it and squeezed firmly down. If this operation be 
carefully performed it is possible to fix the wires in their positions, though 
the thin platinum wires gave some trouble owing to their great fragility. 

D 2 
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When the apparatus is used for testing the conductivity of powders, the 
lower shell is first filled and the wires placed in position on the upper 
surface of the substance, the current turned on and the heating wire drawn 
tight. The upper shell is now laid above the first, and filled through a 
slot made along the top, the powder being distributed to fill it completely. 

The apparatus is shown generally in fig. 1. Fig. 2 is a horizontal section 


Fig. 1. 




in the plane of junction of the two half cylinders; A A, BB are the two 
thermometric wires, and CD the heating wire fixed to brass plates. The 
length of the cylinders used was 15 cm., and the diameter 9 cm. The 
heating wire had a length of 25 cm. in the later experiments. 

3. Thermometric Arrangements .—(1) When the difference of temperatures 
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1904.] Finding the Conductivity for Heat 37 

was measured by the change of resistance of platinum, a piece of platinum 
wire was coated with copper and the whole drawn out till the platinum had 
a thickness of about 0*035 mm., that of the copper being about 0*27 mm. 
Any change of resistance was thus appreciably due to the platinum. Two 
lengths from the same piece were taken, laid together, and being connected 
with soft wax, were bared as nearly as possible to the same extent of about 
3 to 4 cm., so that their resistances were very nearly equal. 

If R 0 be the resistance of the inner one, R' 0 that of the outer before 
heating, let 

E 0 : R' 0 = a: 

when compared by the Wheatstone’s Bridge method, taking, say, a = 1000 
and /3 a number not far from 1000. After heating, and when the 
temperature has become steady, 

R 0 becomes R x = E 0 [1+c —0 O )], 

R' 0 „ R , i= R r 0 [1+c (0 2 —0 O )], 

where c is the temperature coefficient of the resistance of platinum, and 0 Q 
is the temperature of the room. 

Suppose now that R x : R^ = a : Then 

l±£jgiZLgo) = £L 
1 + c (0 2 —0 O ) ft l 

If the total changes of resistance are only comparatively small fractions of 
the original resistance, we have 

HOi-e,) = 08-/SO/A. 

Of course, if c (0 3 — 0 0 ) is a considerable fraction, the resistances must be 
found separately. In the experiments actually made, O 2 —0 o amounted to 
about 30° at most, and c = 1/400, so that the assumption here made is 
sufficiently accurate. 

(2) So much difficulty was experienced in the use of these fine platinum 
wires, especially in testing loose powders such as sand or mould, that they 
were ultimately given up, and the difference of temperatures found by the 
thermo-electric junctions already described. The arrangement used was a 
wire of German silver of 0*2 mm. diameter between two iron wires of the 
same diameter, stretched parallel to the axis of the cylinder so that the two 
wires, and, therefore, the junctions, were at the distances apart previously 
chosen, that is to say, 1 cm. and 3 cm. from the heating wire. As the 
resistance of the galvanometer employed was about 800 ohms, the deflection 
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was sensibly proportional to the difference of temperature. By independent 
experiments one scale division represented 0*70° C. 

4. With regard to the choice of the distances a , b from the central wire, it 
is clear that these should not be taken too nearly equal; for, unless & be a 
moderate multiple of a , any error in the measurement of one of them, say a, 
will introduce a large error in log,, bja. 

Let a f for instance = a (1 — 1/%), 

log* bja' = log b/a + l/n , nearly. 

Thus, if a = 1 cm., b = 1*5 cm., and - a == 0*5 mm., 

n 

log b/a = 0-405, l/% = 0*05, 

so that the error is about 12 per cent, of the whole. 

By taking = 3 the error lies between 4 and 5 per cent. 

As between the two pairs of distances 1, 6 and 2, 6 cm., we should rather 
choose the latter, as log* 6 is not quite 2 log e 3, and an error of 0*5 mm. in 
2 cm. is only half as important as the same error in 1 cm. 

In the present form of the apparatus I have therefore taken a — 1 cm., 
b = 3 cm.; and in a larger model, a = 2 cm., b = 6 cm. 

5. It is desirable to know how long an experiment need be continued to 
get a fair approximation to the steady state. For this purpose it is necessary 
to know how the temperature rises at any point. 

In the second part of this paper the solution is given of the problem of the 
diffusion of heat in an infinite cylinder from a steady source at the axis. 
Taking the result as a guide for the present case, we have 

n rbj'Snt 

= f U -'e-^clx, 

" J a/ */nt 

where x 3 == r 3 /%£, and r is the distance from the axis and n the diffusivity. 

Assuming, provisionally, that r 3 /4%£ is large enough to make r 3 /4%£ small, 
we may expand e~~ r ^ hlt , and retain only the first two terms. Thus 


The first term of the expression gives the temperature of the steady state, 
and if the temperature be 5 per cent, below its final value, 


& 3 —• a? 
8nt 


nr log, - 




If a = 1 cm., b = 3 cm., 


t = 18*2/%. 
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In the case of sand, for which n is given, 00032 (Everett's ‘Units and 
Physical Constants ') 

t = 5702, or about 1^ hours. 

To get within 1 per cent., we should have to continue the heating for about 
7 J hours. 

It may, perhaps, be worth remarking that the diffusivity may be found 
directly from the expression for ddjdt given in the second part of the paper 
and quoted above ; for 

dt t 

which is a maximum when 

t = r^/4n. 

Taking r = 3 cm., the time required is only about one-eighth of that 
previously necessary. 

6. As an illustration of the results which the method gives, I add a list of 
some substances whose conductivities have been found by it. In the case of 
the woods, the flow of heat is across the fibres. In the second column the 
letter E indicates that the difference of temperatures was found by the change 
of resistance of the platinum wires, J that it was found by thermo-electric 
junctions. The distances of the thermometric wires from the central wires 
are given in the third and fourth columns; in the earlier experiments, bja 
was taken nearly = 1*5, while in the later ones bja is about 3. The 
numbers agree fairly well, though, for the reasons given above, I think the 
latter ratio preferable. The fifth column shows the size of the platinoid wire 
used as the source of heat, and the sixth the current passing through it. 

Some determinations by other observers are added; these are taken from 
Everett's ‘ Units and Physical Constants.' 

I wish to express my indebtedness to Mr. William Mitchell, M.A„ by whom 
many of these determinations were made, for this and other assistance he has 
given me in the course of the work. 
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Authority. 

G. Forbes. 

Peclet. 

B.A. Committee. 

(Ordinary) Peclet. 

(Fine) Peclet. 

G. Forbes. 

B.A. Committee. 

Peclet. 

G. Forbes. 

(Mahogany) Peclet. 

Conduc¬ 
tivity 
x 10 4 . 

^ »0 CO NCOOOOt- 04 00 O i-h O CO 

OOJO 5 OOCDO 50 O»OOi>i>O(M(M>OOHCOCOlOCOb-(MOO 

(MCOiMGXMOlCOCO^COCOOiMOKMOCDJxHOJt-^HHH 
r-H H t-H 

Currents 

in 

amperes. 

CO 1 — CD H (M t- (M 05 N kO CO 

000005 00 CO <M 00 (M 00 00 GO Cl O 04 CM | | 

O O O O O T-H o T-H o O O r-H r-H r-H r-H 

No. of 
wire 
S.W.G. 

^ O O r^O^O^OO ^ I 1 

CO hH CO 'rH CO ^ CO ^ 'rH CO CO COCO | | 

Distances from the 
central wire. 

rO 

Jt— J>» CO r-H CO N 

10^0 »0 IQ 05 lO 05 CO ^ CO 05 GO 05 00 J | 

r-H r-H CO r-H r-H 04 r-H 04 t-H r-H 04 04 04 04 04 


05 05 lO i—CD C 4 Jt>- 

r-H 05 05 000000500 04 05 O 05 05 | | 

r-H O O r-H r-H r-H r-H r-H O O r-H O r-H OO 

Thermo¬ 

metric 

means. 

P3P3P3 P3 P5 H P^ 1-5 p? P3 P3HH I I 

Substance. 

Fir. 

Butternut . 

Mahogany .. 

Beech ... 

Plaster-of-Paris (made-up) 

Sand (fine). 

Garden mould (dry). 

Sawdust (fir). 

j> . . 

» . 
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Part II. 

In connection with the foregoing method, based on the steady state, it seemed 
desirable to know the mode in which this state is attained at any time after the 
heating is begun. As a preliminary, the problem of the heating of a rod by 
heat supplied uniformly in the middle, was first attacked in Fourier’s way by 
treating it as the limit of a ring of infinite radius. The results, though 
complicated in form, may be worth recording, as the conditions approximate 
to those actually met with when one end of a metal rod is put into a hot 
flame. 

The case of a small spherical cavity is added, and may be taken as an 
illustration of the way in which the heat given out by a small morsel of radium 
is propagated through an infinite solid. 

In all the questions treated in this part the body is supposed infinite in the 
directions in which the flow of heat is considered, the ring and thin spherical 
shell excepted. 

Case 1. Thin Circular Ring with Surface Impermeable to Heat .—Consider 
first the case of a uniform thin circular ring heated at one spot. 


Let k = conductivity, <r = specific heat per unit volume, 
A = cross-section of ring, c = its mean radius, 

s = an arc of ring measured from the heated spot, 

= the angle subtended by s at the centre, so that s = c<£, 
H = heat supplied at the spot per second. 

The general equation of conduction is 


7 a a &0 


or 


k<P0 d£ 
c 3 defy 2 a dt 


in 


6 being the temperature. 

If the heat has been supplied for a long time t, the temperature will be 
sensibly uniform throughout and equal to 


Hi / 2itcA(t. 

If we include the effect of conduction we may satisfy Equation (1) by adding 
to this a term of the form B$ 3 -f C$-f D, where 

B = -Re^irAk 

and the other constants are undetermined. 

To complete the general solution we have to add terms of the form 


L + SMfCOs i<j>e~ xt .... (2), 

where \<r = ki^Jc 2 .. (3). 
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We have also to satisfy the conditions that the flow of heat outwards at the 
origin 0 is constantly equal to |H on either side, that is, corresponding to 
<f> = 0, <£ = 2i r. 

All these conditions are satisfied by the expression 

6 - ;r S!_+-^ 7 ((/.-7r) 2 + L + SM i cos^fi- w . (4). 

2ttcAct 27tA k^ ' w 

It will be noted that the last terms give no flow at 0, on either side, nor at 
the diametrically opposite point, as ought, from symmetry to be the case. 

The constants L, M are determined in the usual way by taking 0 = 0, when 
t = 0, all over the ring. Thus 


And, finally, 


J o J o 

r 2 ir TT r2iT 

Mi j cos 2 i<f> defy + — j 7r) 3 cos i<j> defy = 0. 

Jq ' Jo 


e = _3L | He 7 r) a _ 7tHc He - 1 cog • , g _ AS 
2ircA(r 4 zttAJc 12A& ttAIc i i 3 

Case 2. Infinite Straight Rod Heated at the Middle. —By putting 
ccf> = x, ijc = 3, 1/c = dz, and c = oo , 

we pass, in Fourier’s way, to the case of an infinite straight rod. 

0 = — — -5- f i cos zx e- w l°dx. 

2irAk 67tA k ttA k J z l 
J o 

The second term, which is infinite, may be treated thus— 
t r 3 /6 = l/l 3 +l/2 3 + l/3 3 + . . . 

~ 2 dz. 


(5). 


C7T 3 /6 = Isc 3 /^ 3 = f z 
J n 


If we integrate the last term by parts, dropping the factor H/7 tA& and 
include this term, we have 

x sin zx . 2kt{ 

■ z . 

J o 


B 


^ -f icos 2# 0- 


+ 1 


6- k *t/<r ( fa + 


0 


2Mf 

* Jn 


cos zx e~ hiz2 t <T dz ) 


and the part outside the integral sign now vanishes at both limits. 
If we put n = k/a, 


cos zxe^ ,lt *dz 


2\/nt 
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f sma e _ al! , & 



and thus finally 



( 6 ). 


Case 3. Bing with Badiating Surface .—If P be the perimeter of a cross 
section, and h the constant of radiation, the equation of the flow of heat along 
the ring is 


If the flow of heat from the source were continued for a long time t, the 
final temperature, neglecting conduction, would be given by 


Including conduction, it is easy to see that the general Equation (7) is 
satisfied by an expression of the form 

0— — e -KPt/A v + I; (*-*) + g-m«,-,r)) + L + COS i(f> e~ M ... ( 8 ), 



where 


(9). 


The condition that the flow at either side of the origin should be equal to 
H/2 requires that 


—B. m Jc. IcA (e~ r,m —e vm ) = H/2, 


when t — 0, 0 = 0 everywhere. 


Thus 



+ B (e m «’-*>+ e~”> <*-">) + 2M; cos i(f> — 0. 


Determining the constants L, M ; in the usual way. 



Thus, finally. 


0 


2-nchV 2 m&k e m "—e- mn 


H MWA, ■ H c JL_ *■<»-) + «—»-) c _ hMAv 

. 7. A sJYYlrr ^ _rwrr 


g—JiVtfA<r _|_ 


„-AP«/Acr He £ COS i<f> _jj/„ »2/c» ( 

•n-Aifc i a + wi 3 


(10). 
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Case 4. Infinite Radiating Rod .— To pass to this case, put c<f> = x, c — oo. 


e m (</> - 7T) _|_ e -m ((f) tt) e -px _|_ e vx e -2cpir 


If then^ 2 = PA/A&, m = cp , and 


= e-w, 




(ii). 


This integral may be expressed in terms of two integrals of the error- 
function type, but it is unnecessary to give the reduction here. 

Case 5. Thin Spherical Shell .—The case of an infinite cylinder heated from 
a wire passing along its axis is the same as the two-dimensional problem of 
a plate heated at a point. 

To treat this, consider first the case of a thin spherical shell, heated steadily 
from a small circular hole, and take the boundary of the shell to be impermeable 
to heat. 

The general equation of the flow along the shell is 




a being the radius of the shell, p = cos S, & = angular distance from hole. 
If e be its thickness, the final state, after a long time t , is given by 


TT/ 



To satisfy Equation (12) generally, including this term, we require a term q 
such that 



whence (.1— p 1 ) ^ ~ - 



And the complete value of 6 will be 


00 


6f+q + t 


0 


where A == i(i + l)&/a 2 cr, P* being the usual zonal harmonic of degree i . 
The flow across any small circle S is 



dp 


The part of this depending on Pi vanishes when S = 0 or p =. 1; and, if 
C == 1, the part derived from q will be equal to H. Thus 


dq 


a 
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The complete expression for 6 is therefore 

To determine M f we have, when t = 0, 


Now 

and 

Hence 


f +i r+i 

Pi 2 cZ/i + j gVid/M = 0. 


d_ 

dfi 


The first term on the left hand vanishes except when i = 0. On the 
hand (1— f) Pj = a (1 + ft) Pi = 2 a, between the limits; and 

(flit 

(1-V) q = (1-/0 log (1-/0 x(l+<0 

which vanishes both when //. = 1 and when n = — 1. Thus 

Mi= -o(2t+l)/*(t + l). 

The constant term M 0 is given from 


M„ 


*oJ dp. — «|log(l— p)dp = —«{(1—/01og(l—/0~(1—/0}tl 

= a (2 log 2—2), 
or M 0 = a (log 2—1). 

Thus the complete expression for 6 is 

m 


6 = 


where 


-—logf+1T —«2Pie~ xs ...... 

4:Tra?e<T l \ 2 / / 

a — H / Aireh. 


When t = 0, 0 = 0 everywhere, and this requires that 

l og Hp +1 + S <f+I) P,= ’°. 

an expression which may be readily verified from the expression 

(1—2 fix -f- == 1 + SP 
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Case 6. —Infinite Plane Plate. This case, which also includes that of an 
infinite cylinder heated along its axis, may be deduced from the foregoing 
one by putting 

l-fi = 2 

r being the distance of any point from the centre, and by taking 
i/a = z, 1/a = dz, X = kz^/cr ; 
also Pj = J, 0 (sr), 



and the expression for 0 is 


0 


_ _ Hj_ r . ,V . H- 


Azrk \ 


21q«s+1 + 


Sh«-f 


z (rz) e~~ nts *dz 


(15), 


where H x stands for the heat supplied per unit length in unit time, and 


n = kjcr. 

The integral, which itself is infinite, is coupled with an infinite term, the 
difference between the two being the finite part on which the solution 
depends. The same point came up in the deduction of the solution for the 
infinite rod from that of the ring, but in that case it was possible to assign a 
finite form for the difference. 

The form of the integral in (15) suggests, as will be seen later, another way 
of obtaining the solution, but it is possible to deduce from it a form of 
solution satisfying all the conditions of the problem. For if we differentiate 
the above expression with regard to the time, we have 


cW _ 
dt 27 rk 


n 


zJ 0 (rz) e~ ntz2 dz ; 

* n 


or, putting 
we have 


f = z\/ nt , x = r/\/ nt , 


d6 = H x 1 

dt 47 rk t 


u, 


u 




It is easy to verify by direct differentiation that 


and therefore 


du . n 
z — +ux = 0, 

dx 

u — Ae- x ^ 


To find the constant A, we observe that when 

x = 0, u = f = |; so that A = \ 

J o 

Thus ^ 0 = (rh-^hlt. 

dt Anrk t 4:irk | 
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But 


a? = r^/nt, and therefore 2dxlos J tdtjt = 0. 


Thus 



(16). 


Another Solution of this Case. —The expression just given for 0 was 
originally obtained by a different method, which may be now given. The 
integral contained in the Equation (15), viz., 

f z~ l J 0 (rz ) e~ %u2 dz, 

may be written, using the same symbol £ as above, in the form 


'Vjo (fr/v'nO 


0 


which shows that 0 is of the form cfr^/nt). 

But the general equation of the flow of heat in two dimensions is 


m id0 

dec 2 dy 2 ndt 3 


which may be satisfied by 0 = <f>(v), where v — (f + y^/nt, if 



v ^ = Ae ~ v/4 . 


and therefore 


The total flow across the boundary of the circle of radius, r, in 


= — 2-rrrk ^ = — 4:7rJcv ^, 
dr dv 


and when r = 0, or v = 0, this is to be equal to H 1 . Therefore 

A == -Hf4irh 



Thus 


the upper limit of integration being taken to make 0 everywhere = 0, when 


t = 0 


Writing for v, this expression becomes the same as that previously found 
in Equation (16). 
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Case 7. Case of an Infinite Solid with small Spherical Source .—The equation 
of motion in this case is 

dd 


h 


f6 + 2cW 

clr 2 r dr 


dt ' 


or 


dHff) _d(Or) 


dr 2 dt 

where n = h/a. 

When t is infinite and the flow has become steady, 

6 = A/r + B, where IcAttA = H. 
But the general Equation (17) is satisfied also by 


(17), 


6r = 


pGO 00 

= <f) (z) cos zr e~~ ntz *dz, and by j <f> (z) si 

Jo 


sin zr e~~ nU2 dz. 


If we choose the latter, and put $(%) = C/z, we may take as the solution 
applicable to the present case 


H 1 , C fsin 


0 = —-f- 

4:7r/c r r 


" A 


t zr 


e~ ntz *dz. 


For, 


(2) When t = 0, 0 = -A +C 
v ' r\4:irk 


(1) When t = oo, it reproduces the steady state; 

f-^4 
= ;(is + f c )' 

C = -2/7r.H/47r*; 


which vanishes if 
(3) -4 tt= H—— f"S3EL*’ e - nt *dz+ —Ceos zre~ nt *dz, 

dr 77* I Z 7T J 

J 0 0 

which is equal to H when r = 0. All the conditions are therefore satisfied. 
The solution therefore is 


0 = JL ifl-i p!HL2T e- nt *dz) = JL l(l — j=.\ r e 
vrrh j* \ 7r J z / 47 tk r ' \/nt J 



















